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ABSTRACT 

Wc present high resohition i?-band images of the central regions of 67 early-type 
galaxies obtained with the Wide Field and Planetary Camera 2 (WFPC2) aboard the 
Hubble Space Telescope (HST). This homogeneously selected sample roughly doubles 
the number of early-type galaxies that have now been imaged at HST resolution, 
and complements similar data on the central regions of radio galaxies and the bulges 
of spiral galaxies. Our sample strikingly confirms the complex morphologies of the 
central regions of early-type galaxies which have become apparent from previous 
studies with HST. In particular, we detect dust, either in the form of nuclear disks 
or with a filamentary distribution, in 43 percent of all galaxies, in good agreement 
with previous estimates. In addition, we find evidence for embedded stellar disks in a 
remarkably large fraction of 51 percent. In 14 of those galaxies the disk- like structures 
are misaligned with the main galaxy, suggesting that they correspond to stellar bars 
in SO galaxies. We analyze the luminosity profiles of the galaxies in our sample, and 
classify galaxies according to their central cusp slope. To a large extent we confirm the 
results from previous HST surveys in that early-type galaxies reveal a clear dichotomy: 
the bright ellipticals {Mb ^ —20.5) are generally boxy and have luminosity profiles 
that break from steep outer power-laws to shallow inner cusps (referred to as 'core' 
galaxies). The fainter ellipticals, on the other hand, typically have disky isophotes and 
luminosity profiles that lack a clear break and have a steep central cusp (referred to 
as 'power-law' galaxies). The advantages and shortcomings of classification schemes 
utilizing the extrapolated central cusp slope 7 are discussed, and it is shown that 
7 might be an inadequate representation for galaxies whose luminosity profile slope 
changes smoothly with radius rather than resembling a broken power-law. Thus we 
introduce a new, alternative parameter, and show how this affects the classification. In 
fact, we find evidence for an 'intermediate' class of galaxies, that cannot unambiguously 
be classified as either core or power-law galaxies, and which have central cusp slopes 
and absolute magnitudes intermediate between those of core and power-law galaxies. 
It is unclear at present, however, whether these galaxies make up a physically distinct 
class or whether distance and/or resolution effects cause them to loose their distinct 
core or power-law characteristics. 

Subject headings: galaxies: elliptical and lenticular, cD — galaxies: nuclei — galaxies: 
structure. 
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1. Introduction 

Ever since the seminal study by Davies et al. (1983) it is known that ehiptical galaxies can 
be divided in two distinct classes. The bright ellipticals {Mb ^ —20.5) generally have a boxy 
appearance, are pressure supported, and often show some form of radio and/or X-ray activity. Low 
luminosity ellipticals, on the other hand, often have disky isophotcs, arc rotationally supported, 
and generally lack significant amounts of activity (Bender 1988; Bender et al. 1989; Nieto et 
al. 1988). Before the Hubble Space Telescope (HST) came online, very little was known about the 
nuclear properties of these systems. Based on ground-based observations, with a typical seeing of 
^ 1 arcsec, ellipticals were considered to have constant density cores with a size that correlates 
with global galaxy properties (e.g., Lauer 1985; Kormendy 1985). High-resolution imaging surveys 
with the HST, however, revealed that the central brightness profiles of early-type galaxies are 
well approximated by a power-law, I{r) oc r~'^ with 7 > (e.g., Lauer et al. 1991, 1992; Crane 
et al. 1993; Ferrarcsc ct al. 1994; Kormendy et al. 1994; Lauer et al. 1995; Carollo et al. 1997a). 
Furthermore, the distribution of 7 was found to be bimodal and to follow the dichotomy of global 
properties inferred from ground-based observations (Jaffe et al. 1994; Lauer et al. 1995; Faber et 
al. 1997) and to be correlated with the global X-ray emission (Pellegrini 1999). Bright ellipticals 
have luminosity profiles that are well fit by a double power-law with a break radius of typically a 
few hundred parsecs. The inner power-law slope has values of < 7 ^ 0.3. This contrasts strongly 
with the low-luminosity ellipticals whose luminosity profiles lack a clear break radius and have 
central power-law profiles with 7 ^ 0.8 on average (Jaff'e et al. 1994; Ferrarese et al. 1994; Lauer 
et al. 1995; Byun et al. 1996; Gebhardt et al. 1996; Faber et al. 1997; Quillen, Bower & Stritzinger 
2000). Following the nomenclature of Lauer et al. (1995) we refer to the former as core galaxies 
and to the latter as power-law galaxies. 

In addition to the central density cusps, the images from the HST have revealed that a 
large fraction of carly-typc galaxies harbor significant amounts of dust (~ 10'^ — lO'' M©), either 
in the form of a nuclear disk, or of a more complex, filamentary or patchy morphology (Jaffe et 
al. 1994; Lauer et al. 1995; van Dokkum &: Franx 1995). In several cases, nuclear stellar disks with 
scale-lengths as small as ~ 20 pc, are found (e.g., van den Bosch et al. 1994; van den Bosch, Jaffe 
& van der Marel 1998; Scorza & van den Bosch 1998). Both the gaseous and stellar disks have 
been shown to be powerful tools for determining the central densities of their parent galaxies, and 
for detecting massive black holes (Harms et al. 1994; Ferrarese, Ford Sz Jaffe 1996; van den Bosch 
& de Zeeuw 1996; Macchetto et al. 1997; Bower et al. 1998; Cretton & van den Bosch 1999; van 
der Marel & van den Bosch 1999). Mainly because of these nuclear disk structures the number of 
galaxies with solid detections of massive black holes (BHs) has increased strongly over the past 
years, and we have reached the point where we can start to address the demographics of BHs (see 
Kormendy & Richstone 1995, Ford et al. 1998, Ho 1999, van der Marel 1999a, and Macchetto 1999 
for recent reviews). 

Black holes also play an important role in connection with the observed central cusps. Two 

different scenarios have been proposed that rely on the presence of massive black holes to explain 
the correlation between cusp slope and luminosity. Faber ct al. (1997) have suggested that 
luminous galaxies are the result of mergers, whereby BHs in the progenitors sink to the center and 
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excavate a core (see also Begelman, Blandford &: Rees 1980; Makino & Ebisuzaki 1996; Quinlan 
1996; Quinlan &: Hcrnquist 1997). Van der Marel (1999b), on the other hand, starts with the 
assumption of constant density cores, and shows that the adiabatic growth of central BHs can 
create cusps whose slopes correlate with luminosity as observed. Thus, whereas the origin of 
the relation between cusp slope, luminosity and black hole mass is still unknown, it has become 
apparent that massive black holes are omnipresent in galactic nuclei. 

Despite the enormous wealth of new data, we are still far from a proper understanding of 

the structure and formation of the central regions of early-type galaxies. Not only have the new 
HST observations of these systems changed our view of their central regions, they have also raised 
several new problems: 

• Is the bimodality in cusp slopes real? Data on elliptical galaxies with kinematically distinct 
cores seems to suggest a much more continuous transition (CaroUo et al. 1997a). 

• What is the origin of the correlation between cusp slope and global characteristics of 
early-type galaxies? 

• What fraction of early- type galaxies contain nuclear disks of stars or dust/gas, and how are 
these structures related to other properties of their parent galaxies? 

• Is there a relation between the central structure of galaxies and the presence or absence of 
nuclear activity? 

• How does the mass of a central BH correlate with other properties as bulge mass and cusp 
slope, and what role do BHs play in determining the properties of the central regions of 
galaxies? 

In order to address these questions one needs a large and unbiased sample of early-type 
galaxies that have been imaged at comparable resolution. On the order of 80 early-type galaxies 
have been imaged with the HST, most of them before the refurbishment mission. We have 
obtained HST images with WFPC2 of the central regions of 67 early-type galaxies as part of our 
HST snap shot program ^ 6357. This roughly doubles the total number of early-type galaxies that 
have been imaged with the HST. Together with the data taken so far and with similar samples of 
spiral bulges (Phillips et al. 1996; Carollo et al. 1997b; Carollo, Stiavelli & Mack 1998) and radio 
galaxies (de Koff et al. 1996; McCarthy et al. 1997; Martel et al. 1999; Verdoes Kleijn et al. 1999) 
these images provide a data base for the investigation of correlations of properties such as dust, 
metallicity, colors, activity, and nuclear structure. 

In this paper, the first in a series, we present the data and describe the luminosity profiles and 
isophotal structures. More detailed analyses of the dust, stellar disks, and the central parameter 
relations are deferred to future papers. 
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2. The Sample 

The sample is compiled from the Lyon/Meudon Extragalactic Database (LEDA)^ by selecting 
all early- type galaxies with radial velocities less than 3400 kms~^, an absolute l^-band magnitude 
less than —18.5, and absolute galactic latitude exceeding 20 degrees (to minimize the effects 
of galactic extinction). In total 130 galaxies made it into the sample, of which 68 have been 
observed successfully by HST. Thus, although our sample is by no means complete, it has been 
homogeneously selected. One galaxy, UGC 5467, turned out to be a spiral galaxy and will be 
discarded in what follows. Note that, although we have attempted to discard obvious duplicates 
(WFPC2) from our sample, several of our objects have been observed in other HST programs 
since the time our sample was selected. In most cases, however, a different filter was used than for 
the data presented here. 

Global properties of the 67 galaxies in our sample, taken from the LEDA (see Paturel et 
al. 1997), are listed in Table ||. Absolute magnitudes (column 2) are based on the distances 
listed in column 4, and based on the Virgo-centric infall corrected heliocentric velocities assuming 
Ho = 80kms^^Mpc^^ (see Paturel et al. 1997 for details). In addition, we list the 1.4 GHz 
radio fluxes, taken from the NRAO VLA Sky Survey (NVSS; Cordon et al. 1998), as well as a 
far-infrared magnitudes defined as 

mpiR = -2.5 log(2.58/6o + /loo) + 14.75 (1) 

(cf. de Vaucouleurs et al. 1991). Here /eo and /loo are the q60 ^m and 100 fim IRAS fluxes 
measured by one of us, adopting a similar approach to that taken by Knapp et al. (1989). 



3. Observations and Data Reduction 

3.1. Observing strategy 

We used the WFPC2 to obtain broad band images of the galaxies in our sample, using the 
F702W filter. This filter, centered on A = 6997A and with a FWHM of ~ 148lA, was chosen to 
compromise between maximizing the number of photons and minimizing the effects of dust on the 
observed morphologies. The nuclei of the galaxies were centered on the Planetary Camera CCD 
(PC), which consists of 800 x 800 pixels of 0.046" x 0.046" each. All exposures were taken with 
the telescope guiding in fine lock, yielding a RMS telescope jitter of ~ 3 milli-arcseconds. More 
detailed information on the WFPC2 can be found in Biretta et al. (1996). 

Previous HST/WFPC2 observations of ellipticals have shown that the power-law galaxies 
have typically a central surface brightness in the y-band of fiQ ~ 14.5, whereas core-galaxies 
typically have /io ~ 16.5. Our aim is to obtain a signal-to-noise (S/N) of 20 or larger in the 
center of the galaxies. This prompted us to use a total exposure time of 1000 seconds with the 
analogue-to-digital gain set to 13.99 electrons/DN (DN is the number of counts). With this 



^www-obs. univ-lyonl.fr 
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setting the read-out noise is 7.02 electrons. To facilitate cosmic ray removal, and to guard against 
saturation in the presence of a bright nuclear point source, each exposure was split in two separate 
exposures of 300 and 700 seconds respectively. 



3.2. Data reduction 

The images are calibrated by the HST calibration 'pipeline' maintained by the Space Telescope 
Science Institute. The standard reduction steps include bias subtraction, dark current subtraction 
and flat-fielding, and are described in detail by Holtzman et al. (1995a). 

Subsequent reduction is done using standard IRAF tasks. Using the data quality files and 
the WFIXUP task, bad pixels are corrected by means of a linear one-dimensional interpolation. 
The alignment between the 300 and 700 exposures is determined by comparing the isophote 
centers from initial analyses of the separate exposures. If the misalignment is significant the 300 
sec exposure is shifted appropriately, using linear interpolation. A difference-image of the 700 sec 
and (shifted) 300 sec images is used to check the accuracy of the alignment. Next, the two images 
are combined with simultaneous removal of cosmic rays, using the IRAF task IMCOMBINE. The 
cosmic ray removal and bad pixel correction is checked by inspection of the residual images between 
the cleaned, combined image and each of the two original frames. Sky levels are determined from 
empty areas in the WF images, and accordingly subtracted. The sky-level is found to be small, 
typically ~ 1 DN. The count rates in the broad-band F702W images are calibrated to magnitudes 
in the Landoldt i?-band, as described in Holtzman et al. (1995b), and assuming a V-R color of 
0.68. The magnitudes thus derived are checked using the inverse sensitivity and the zero point of 
the magnitude scale provided by the HST pipeline, whereby it is assumed that ellipticals have the 
spectrum of a K-III giant. The consistency is found to be better than 0.01 dex. 

The 300 and 700 sec exposures are reduced separately as well. Cosmic rays are removed using 
the IRAF task COSMICRAY. These images are used to check the consistency of structures, for 
obtaining estimates of the uncertainties due to random noise, and for cases where saturation has 
affected the 700 sec exposure (see § ^ ). 

Despite the high spatial resolution of the refurbished HST, the central few tenths of an arcsec 
of the galaxies are influenced by the effects of point-spread function (PSF) smearing. Because of 
the stability of the HST PSF, these effects can be corrected for by means of deconvolution with 



the Richardson-Lucy (RL) algorithm (Lucy 1974). Based on tests with model galaxies (see § 4.1 
below) we deconvolve the reduced 300 sec, 700 sec and combined images of each galaxy with 40 
iterations of the RL algorithm. For three galaxies with faint and shallow surface brightness profiles 
(NGC 3613, NGC 4168 and NGC 4365), the RL-algorithm with 40 iterations ends up amphfying 
the noise too strongly for an accurate analysis of the central surface brightness profile, and only 
15 iterations are used in these cases. For each galaxy we use a separate model PSF (5.7" x 5.7"), 
created with the TINYTIM software package (Krist 1992), centered on the location of the galaxy's 
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center in the PC CCD^ Contour plots of the deconvolved images are presented in the Appendix. 



4. Data Analysis 

The main goal of this paper is to analyze the isophotal structure and luminosity profiles of 
the galaxies. For the isophotal analysis we use the isophote fitting task ELLIPSE in IRAF. For 
each isophote, the center, ellipticity (e), and position angle (6) are computed. In addition, the 
third and fourth order deviations of the isophote from a pure ellipse are determined. These are 
described by the amplitudes (a„ and 5„ with n = 3,4) of the sin and cos 3^ and 46 terms of the 
following Fourier expansion (see Jedrzejewski 1987): 

= [an sin {ne) + bn cosine)] (2) 

For a perfectly elliptical isophote these coefficients are all equal to zero. Typically the absolute 
values of a„ and 6„ are small, rarely exceeding 0.02. The 64 coefficient is of special interest, since 
positive (negative) 64 values correspond to disky (boxy) isophotes. 

In order to parameterize the luminosity profiles they are fitted with a double power law, 
known as the "Nuker"-law profile (Lauer et al. 1995): 



Ur) = If, 2(^-^)/" ( - 



'1 
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which has I{r) oc at r" ^ rj^, and I{r) oc r~'^ at r" <C rj^. The parameter a describes the 
sharpness of the break from the outer to the inner power law. The break-radius, r;,, is the radius 
of maximum curvature in log-log coordinates, and lb is the surface brightness at r = rf,. The 
Nuker law has been shown to accurately fit the central luminosity profiles of early-type galaxies 
(e.g., Byun et al. 1996). Note that equation (^) is only intended to fit the luminosity profiles 
in the inner ~ 20" of early-type galaxies. For typical fitted values of (3 there must be a further 
downturn in the profile at larger radii for the total luminosity to be finite. 



4.1. Model galaxies 

In order to test the accuracy of our reduction and analysis we constructed a large set of model 
galaxies. Each model galaxy has a luminosity profile of the form of equation (^), has constant 
ellipticity and position angle, and has perfectly elliptical isophotes (i.e., a„ = 6„ = for n > 3). A 
fake image is created by integrating the model surface brightness distribution over each pixel, and 
by adding photon and read-out noise. Pixel sizes are taken the same as for the true data. The 
frame size of 200 x 200 pixels is smaller than the actual HST images, but satisfies the requirements 
for our testing purposes. 



^The exact shape of the HST PSF depends on its location in the focal plane. 
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The resulting image, refered to as 'unconvolved' image, is convolved with the HST PSF (the 
'convolved' image), and subsequently deconvolved using A^iter iterations of the RL algorithm (the 
'deconvolved' image). The PSFs used for the convolution and deconvolution are constructed with 
TINYTIM, and differ slightly from each other in order to mimick errors in the exact PSF shape. 
The unconvolved and deconvolved model images are analyzed in the same way as the data: an 
isophotal analysis is performed, and we fit the Nuker-law profile to the major axis luminosity 
profile determined from the isophote fitting. The deconvolution is tested with 10, 15, 25, 40 and 
80 iterations. The best results were obtained for A^iter = 40, in good agreement with Lauer et 
al. (1998). However, for galaxies with shallow central surface brightness profiles A'^iter = 15 already 
suffices to accurately recover the unconvolved image, and more iterations only amplify the noise. 

In total we have constructed 15 x 50 model galaxies. Each of the 15 sets of model galaxies is 
characterized by a value of 7, ranging from 7 = to 7 = 1.4 in steps of 0.1. For each set 50 model 
galaxies are created whose luminosity profile parameters are varied at the few percent level. The 
ellipticity and orientation of each model galaxy are chosen randomly from the intervals [0, 0.5] 
and [—90,90], respectively. In addition, the center of the galaxy, in fractional pixels, is drawn 
randomly as well. 

In Figure |^ we plot the isophotal parameters of one of the model galaxies with 7 = 0.7. Even 
though this model intrinsically has e = 0.26 over its entire radial range, the ellipticities found by 
the isophote-fitting routine are significantly lower for r ^ 0.2", even for the unconvolved image 
(solid triangles). In addition, the higher-order parameters become noisy and reveal fluctuations 
inside this radius. These effects owe mainly to the discrete pixel sampling and the sub-pixel 
interpolation used by the isophote-fitting routine, and indicate that one can not trust the isophotal 
parameters for r ^ 0.2" (~ 4 pixels). Note, however, that the isophotal parameters of the 
deconvolved model galaxy are remarkably similar to those of the unconvolved model, indicating 
that the deconvolution routine used does not introduce spurious artefacts. 

In order to examine the accuracy with which we can recover the Nuker parameters of the 
luminosity profile used as input for the model, we first focus on the unconvolved images. For each 
model i we compute 

A//i(r) = fio^i{r) - iJ,iso,i{r) (4) 

with /io,i(^) the input surface brightness profile (in magn. arcsec"^) of equation (^), and Hiso,i{r) 
the major axis lumninosity profile determined from the isophotal analysis of the unconvolved 
image. For each of the 15 sets of model galaxies (that have the same value of 7), we then compute 
the average (A^)(r) and the corresponding standard deviation a^^iir), where the average is over 
the 50 models in each set. 

Figure ^ plots (A//)(r) and (T^^{r) for 8 of the 15 sets. For r ^ 0.2" the luminosity profiles 
derived from the isophotal analysis are in excellent agreement with the intrinsic profile. For 
smaller radii, however, the error increases strongly, especially for the profiles with steep central 
cusp slopes. Even more worrysome is that the average residual (A/i)(r) strongly deviates from 
zero at small radii, implying that the main contribution to the error is of systematic rather than 
random nature. This is a reflection of the problems with the isophotal analysis due to the discrete, 
pixelized nature of the data and the sub-pixel interpolation. As Jedrzejewski (1987) already 
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pointed out, the sub-pixel interpolation scheme used can introduce spurious artefacts especially 
when the gradients in the luminosity profiles are large (i.e., when they are steeply cusped). 

Clearly, when determining the intrinsic central cusp slopes of galaxies one can not use the 
surface brightness profiles determined from the isophotal analysis. We therefore use a different 
approach. For r > 0.2" we fit the Nuker law profiles (equation [^) directly to the major axis 
luminosity profile determined from the isophotal fit to the image. For r < 0.2" , however, we 
extract the actual pixel values of the image from the data, which we compare to the model 
predictions given by 

ht{x,y)=JJ I{m)dxdy, (5) 

pixel 

with I{r) the Nuker law profile and 

Here {x, y) are the coordinates centered on the galaxy center {xq, yo) with the x-axis rotated by an 
angle 9q with respect to the pixel coordinate system, and eo and ^0 the ellipticity and position 
angle of the isophotes inside 0.2", respectively. We make the assumption that xq, yo, eo, and are 
all constant inside 0.2" and we compute their values as the average of the corresponding isophotal 
parameters in the radial interval from 0.2" to 0.4". In cases where the isophotes in this radial 
range are not representative of the galaxy, e.g., distorted by dust, a different more appropriate 
radial range is chosen. 

The advantage of this method is that by using the observed pixel values to determine the best 
fit Nuker profile we do not suffer from any problems with sub-pixel interpolation which hampers a 
proper isophotal analysis at small radii. The disadvantage, however, is that if the actual isophotal 
parameters vary with radius inside 0.2", this will introduce systematic errors in the parameters of 
the best fit Nuker law profile. However, detailed tests, described below, show that these errors are 
much smaller than those stemming from sub-pixel interpolation. Note that a somewhat similar 
scheme was used by Byun et al. (1996), who used the total integrated light inside 0.1" as an 
additional constraint for the Nuker-law fitting procedure. 

Figure ^ illustrates the advantage of our new method. The solid lines correspond to the 
intrinsic major axis luminosity profile of the same model as in Figure || (i.e., with 7 = 0.7). The 
upper profiles with the filled circles correspond to the luminosity profiles determined from the 
isophotal analysis. As is evident, inside ~ 0.1" the sub-pixel interpolation used by the isophotal 
analysis introduces systematic errors. Note again that the results for the unconvolved (left 
panels) and deconvolved (right panels) model images are virtually identical, indicating that the 
deconvolution procedure used can accurately correct for the PSF smearing. The lower profiles 
with open symbols, offset by 1.5 magn, correspond to the luminosity profile determined with the 
new method. The open circles (r > 0.2") are the luminosities determined from the isophotal 
analysis. The open squares correspond to the actual pixel values of the model image inside 0.2"g. 



^Since in general the (center of the) pixels do not lie exactly along the major axis, we make use of the following 
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As is evident from the lower panels, this new method allows us to recover the intrinsic luminosity 
profile with high accuracy down to very small radii. 

In order to determine the parameters of the Nuker law profile that best fits the data we 
use a x^-iiiiiiimization technique. For a given set of parameters (a, P, 7, rb, Ih), we determine 

= Xi + xi- Here Xi ^^'^ X2 x^-values of the pixel-value fitting inside 0.2" and the 

luminosity profile fitting at r > 0.2" , respectively. As errors we use Poisson statistics for the pixel 
values and the nominal errors given by the isophotal analysis, respectively. In practice, these latter 
errors are underestimated, as is evident from the fact that the normalized xi of the best fit model 
is typically significantly more than unity. We therefore introduce the parameter /, which scales 
the nominal error given by the isophotal analysis, and which thus sets the relative contributions 
of Xi and xi to the total x^ of the fit. In the left panel of Figure |^ we plot the average difference 
between the intrinsic cusp steepness 7modei and the best-fit Nuker-law profile 7fit, as function of / 
for the set of models with 7 = 0.5. The averages are taken over all 50 model galaxies in this set. 
As expected, the fits are better for / > 1. A value of / ~ 5 yields values for 7fit in good agreement 
with 7modei- For / ^ 1, too much relative weight is given to the pixel values inside 0.2", and 7fit is 
too large. For too large /, x^ is completely dominated by the isophotal luminosities outside 0.2" 
and the central cusp slope becomes poorly constrained, as is evident from the large errorbars. In 
what follows we therefore adhere to / = 5, for which in general we find a normalized of the 
best fit model close to unity. 

We have constructed luminosity profiles for all 750 deconvolved model images using our new 
method, and determined the best fit Nuker law parameters using the x^-minimization technique 
described above with / = 5. In Figure |^ we plot jfn versus 7modei for all 750 models. As is evident, 
this method of analyzing luminosity profiles allows an accurate determination of the central cusp 
slope, and we use the same method to analyze the luminosity profiles of the actual HST images. 



5. Results 

5.1. Isophotal Analysis 

For each galaxy there are six reduced images: the non-deconvolved and deconvolved 300 sec, 
700 sec, and combined images. An isophotal analysis on each of these six images is performed. All 
images are checked for contamination such as cosmic rays that have not been properly removed, 
foreground stars, dust filaments, and bright globular clusters. If necessary, these objects are 
masked out. Ten galaxies are so severely influenced by dust that meaningful estimates of their 
luminosity profiles and isophotal structures can not be obtained (all galaxies with filamentary 
dust of "type" III as well as NGC 4233 and NGC 4494, which have a large dust disk, see § ^ 



way to visualize the fit residuals in the luminosity profile plot: We plot the pixel values at the radius Vpiot defined 
such that I{rpiot) = Ifit{x, y), i.e., the difference between the pixel value and the Nuker-law profile at Vpiot is equal to 
the difference between the actual pixel value and the prediction Ifuix, y) of the best-fit model. This way the residual 
is accurately visualized. 
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below). In most of the following discussion we discard these galaxies and focus on the sample of 
57 galaxies for which the HST images provide a meaningful picture of the structure of the stellar 
component. In what follows we refer to this sub-sample as the unperturbed sample. 

Isophotal parameters are determined down to r = 0.03", and presented in the Appendix. 
Note that, although we plot the isophotal parameters down to r = 0.03", only the parameters 
for r ^ 0.2" are considered for further analysis. For comparison, the isophotal parameters of the 
deconvolved image (open circles) are overplotted with the ones of the non-deconvolved image 
(solid triangles). At larger radii, unaffected by convolution effects, the latter data, which is less 
noisy, nicely complements the deconvolved data. 

Previous high-resolution imaging with the HST has revealed that the central regions of 
early-type galaxies are very complex environments. This is strikingly confirmed with our sample, 
and is clearly reflected in the isophotal parameters: the isophotal shapes can vary strongly with 
radius and are very different from one galaxy to the other. This makes it very hard to categorize 
the central regions of early-type galaxies in a small number of classes. It has become customary 
to classify ellipticals as either disky or boxy, depending on the sign of the 64-parameter. This 
generally works well for the outer parts of ellipticals, where the changes in the isophotal shape are 
only modest. However, the central regions reveal a vastly more complex behavior. A single galaxy 
can have isophotes that change from disky to boxy and back to disky all within the inner 10", even 
when there are no signs of dust. Clearly, a bimodal classification does not suffice to describe the 
bewildering variety of morphologies apparent in the inner ~ 10". Another approach is therefore 
employed: Rather than classifying galaxies as either disky or boxy a simple letter coding, referred 
to as isophotal code, is used to describe the individual isophotal shapes of each galaxy. The letters 
indicate regions with common isophotal shape, starting from the center at r = 0.2" out to r ~ 25", 
where d, b, 0, x and ? refer to disky, boxy, neither disky nor boxy, undetermined due to dust, 
and undetermined due to small surface brightness gradient, respectively. As an example, the 
isophotal code xOd for NGC 2592 indicates that for small radii the isophotal shape could not be 
determined due to dust (in this case a small dust disk, see § ^]^), followed by a region where the 
isophotes are elliptical, and which become disky at larger radii. As a general rule, a galaxy region 
is classified as disky (boxy), if ^4 > -1-0.01 (64 < —0.01) for several consecutive isophotes while the 
other high-order parameters are not significantly different from zero. Exceptions are misaligned 
structures which also result in non-zero 04 values. A more detailed description of these intriguing 



structures is given in § 5.5 



The isophotal codes for each galaxy are listed in column (5) in Table Clearly, this encoding 
is somewhat subjective. Nevertheless, it is well suited to describe the complex and diverse 
isophotal behavior of the galaxies. In addition to the isophotal code, the median ellipticity (e) 
calculated over the radial interval from 1.0" to 10.0" is listed in Column (8) of Table ^. Note 
that these median values may not be very meaningful for galaxies that have rapidly changing 
ellipticities. Large variations may be real or may reflect the presence of dust, and are easily 
identifled by their relatively large errors (A(e) ^ 0.1). 
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5.2. Luminosity Profiles 

Luminosity profiles along both the major and minor axes of the deconvolved images are 



determined and subsequently fitted with the Nuker-law profile as described in detail in § 4.1. In 
some cases a deviation from equation (^) is visible in the inner ~ 0.1" as excess light, and pixels 
inside 0.1" are ignored in the fitting procedure. These (largely unresolved) 'nuclei' are discussed in 



more detail in § 5/7. In the cases where a central dust disk is present, we exclude any data inside 
the major axis radius of the dust disk from the fit. The best fit parameters are listed in Tables § 
and ^ for the major and minor axes, respectively. Plots of the luminosity profiles are presented in 
the Appendix, with the best fit Nuker law profiles overplotted (solid lines). In cases where the 
luminosity profiles do not reveal a clear change in slope over the entire radial range (NGC 4494 
and NGC 3377 along the minor axis) the observed luminosity profile is fitted by a single power 
law: 

I{r) = hr-" (7) 

For these galaxies, the values listed for correspond to Iq in Tables ^ and |^, and no values for a, 
7 and are given. 



5.2.1. Cusps versus Cores 

As discussed in § |^, HST images have revealed that all early-type galaxies have luminosity 
profiles with central cusps. More importantly, it was found that the distribution of cusp slopes is 
bimodal and in support of the dichotomy of early-type galaxies as inferred from their isophotal 
and kinematical structure. In their study of the central parameter relations of early-type galaxies, 
Faber et al. (1997; hereafter F97) classified all galaxies with 7 < 0.3 and r^, > 0.16" as core 
galaxies. All other galaxies were classified as power-law galaxies. It is important to realize that 
such classification scheme suffers from distance effects. Galaxies at large distances will have break 
radii r^, that are smaller than the resolution limit of the simulations. Consequently, the central 
cusp slopes of the observed luminosity profiles of these systems are to be compared to /?, and not 
7, of less distant systems. Similarly, galaxies that are much closer than the main population of 
systems might reveal central luminosity profiles that deviate from the simple double-power law. 
Indeed, F97 pointed out that both M31 and M32 are poorly fit by the Nuker-law profile, and that 
their classification would have been different had they been at the distance of Virgo. 

With these caveats in mind, we now follow F97 and attempt to classify the galaxies in 
our unperturbed sample as either 'core' or 'power-law'. However, when resorting to the same 
classification scheme as used by F97 we are confronted with another shortcoming concerning 
luminosity profiles having small a. Such profiles occur in our sample for two reasons: because of 
the distance-problem alluded to above, or because they intrinsically have luminosity profiles for 
which the slope changes smoothly with radius. Some galaxies in the latter category can still be fit 
by a Nuker-law profile, while for others, even small a Nuker-law profiles fits cannot accomodate all 
profile features. In all of these cases we generally find best-fitting Nuker-law profiles with 7 ~ 0. 
However, 7 is only a proper representation of the gradient in the luminosity profile at radii r for 
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which {r/rt,)" <C 1. For small a, however, {r/ri,)'^ can be close to unity, even when r is significantly 
smaller than rh, and 7 therefore becomes a measure of the gradient of the luminosity profile at 
radii much smaller than the actual resolution limit. Since one has to be sceptical about basing 
a galaxy's classification on extrapolated quantities, we complement our classification scheme by 
introducing the parameter 



7' 



d log 1 1 + 
d log r J, ,=0.1" 1 + 



(8) 



with r' = 0.1". Thus, 7' is determined directly from the best-fit Nuker-law, and is a measure of 
the gradient of the luminosity profile at r = 0.1". The values of 7' along the major and minor axes 
of the luminosity profiles are listed in Tables ^ and ^, respectively. 

In the upper panels of Figure |5| we compare the cusp slopes of the luminosity profiles along 
the major and minor axes. Both 7 and 7' give very similar results along both axes, indicating 
that our method of using the actual observed pixel values inside 0.2" to constrain the Nuker-law 
fit yields a robust measure of the central cusp slope. Comparing 7 to 7' for both the major 
and minor axis luminosity profiles (lower panels of Figure ||) yields a more complicated pixture. 
Galaxies with 7 > 0.5 in general have 7' ~ 7. However, for 7 < 0.3 about half of the galaxies have 
7' < 0.3, and the other half has 7' > 0.3. We therefore use the following classification scheme, 
based on both 7niaj and 7maj0- core galaxies have 7^a^j < 0.3, power-law galaxies have 7maj 

> 0.3, 

and galaxies with 7^3^j > 0.3 and 7maj < 0.3 are provisionally refered to as 'intermediate' galaxies. 
Note that this classification scheme is more conservative in designating galaxies as core galaxies as 
the scheme used by F97: In our sample, only 56 percent of the galaxies classified as core galaxies 
according to the F97 scheme are classified similarly with our definition. The other 44 percent are 
classified as intermediate. Column (2) in Table ^ indicates the classification of our galaxies as 
either core 'n', power-law '\', or intermediate ')'. A question mark indicates that the galaxy is not 
in the unperturbed sample, and that no attempt has been made to classify this galaxy. 

Figure ^ presents histograms of the distributions of 7 and 7'. Although the distribution of 
cusp slopes reveals a clear hint for bimodality, consistent with the results of F97, we find less of 
a deficit of galaxies with 0.3 < 7 ^ 0.5. Most of the galaxies with 7 in this range, and which are 
thus classified as power-law galaxies (e.g., NGC 5576, NGC 5796, and NGC 5831), have small 
values of a and no disky isophotes, similar to the intermediate galaxies. The intermediate galaxies 
themselves follow the same distribution of 7 as the core galaxies, but at the same time occupy 
the I0W-7' end of the power-law galaxies. Thus, whereas we confirm the bimodality in cusp slopes 
found by F97, we find that ~ 15 percent of the galaxies in our unperturbed sample have luminosity 
profiles whose profile slopes change very gradually with radius (reflected by the small value of a), 
making a classification based on central cusp slope somewhat ambiguous (see discussion in § 5.4.2| 
below) . 



6 



our results would be identical had we used and 7n,i„ 
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5.3. Dust Properties 

As discussed in the introduction, the high resolution images obtained with the HST 
have strongly increased the detection rate of dust in early-type galaxies. Two different dust 
morphologies are distinguished: disks (or tori) and filaments (indicated by 'd' and 'f, respectively, 
in Table |2|). The amount of filamentary dust present is indicated by a Roman numeral, ranging 
from I (small traces of dust that hardly influence the isophotal shapes) to III (large amounts of 
dust that prevent a meaningful analysis of the isophotes and luminosity profiles) . Types II and III 
are easily discerned from an inspection of the images by eye, whereas filamentary dust of class I is 
only detected after a more detailed investigation. A detailed discussion of the dust properties in 
our sample is presented in Tran et al. (2000). Here we merely list the main statistics. 

In 10 galaxies (15 percent) we find a nuclear dust disk, with a diameter ranging from 0.3" 
(NGC 4648) to 12.3" (NGC 4233). Most remarkably, 8 of these galaxies are classified as power-law 
galaxies, whereas the two remaining, NGC 4233 and NGC 4494, have a dust disk that is too 
large to be able to classify the galaxy based on its luminosity profile. Filamentary dust is found 
in 19 cases (28 percent). In total, we thus find evidence for dust in 43 percent of our galaxies. 
Van Dokkum & Franx (1995), analyzing 64 early-type galaxies imaged with the HST before the 
refurbishment mission, found dust in 48 percent of the cases and pointed out that a significant 
fraction of the dust will be missed because of inclination effects. They estimated that ~ 75 percent 
of the early type galaxies harbor significant amounts of dust. Our results are in good agreement 
with their estimates. 



5.4. The dichotomy amongst early- type galaxies 

5.4--1. Isophotal shapes 

As mentioned in the introduction, early-type galaxies come in two different classes, core 
and power-law galaxies. The latter are relatively faint {Mb ^ —20.5), typically have disky 
isophotes, and are thought to be rotationally supported spheroidal systems. Core galaxies, on 
the other hand, are relatively bright systems, often with boxy isophotes, that are thought to be 
triaxial and henceforth pressure-supported. We now investigate to what extent the galaxies in our 



sample obey this dichotomy. As discussed in § 5.1, the isophotal structures of the central regions 



are very complex, and we introduced an isophotal code (hereafter IC) to describe the observed 
morphologies. The galaxies can roughly be divided into four classes: disky ('d', ICs with d but no 
6), boxy ('b', ICs with h but no d), disky-boxy ('db', ICs with both d and 5), and regular ('0', ICs 
with no d or b). 

In Table | we correlate this classification based on the isophotal structure with that based 
on the luminosity profiles. We only consider the 57 galaxies in our unperturbed sample. In 
general, our sample nicely confirms the dichotomy that has emerged from previous studies (e.g., 
Jaffe et al. 1994; Lauer et al. 1995): 20 of the 21 disky ('d') galaxies are classified as power-laws, 
while 7 of the 9 core galaxies are classified as either boxy ('b') or regular ('0'). However, not all 
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galaxies fall in either of these two classes. As we indicated in § ^.2.1 , roughly twelve percent of 



our galaxies classify as 'intermediate' based on their luminosity profiles. Remarkably, none of 
these galaxies reveal any sign of diskiness. In fact, two of them classify as boxy, whereas the other 
four are regular. As such, these intermediate galaxies are more similar to core galaxies than to 
power-law galaxies. Of the in total 8 galaxies with both disky and boxy isophotes ('db'), 7 classify 
as power-law galaxies. In total, we thus find that 93 percent of all galaxies with disky isophotes 
are power-law galaxies. If one takes into account that nuclear disks are much easier detectable in 
core galaxies than in power-law galaxies, it is clear that there is a very strong correlation between 
the diskiness and central cusp slope of early-type galaxies. 

Four of the in total 57 galaxies in our unperturbed sample deviate from these general trends: 
NGC 3078 and NGC 5576 are the only power-law galaxies (out of 41) with boxy isophotes. 
NGC 3078 is also the galaxy with the highest 1.4 GHz flux in our sample, and reveals a central 
dust disk with a diameter of ~ 1.4". The dust absorption combined with the presence of a strong 
AGN makes the classification rather uncertain, and it is thus not inconceivable that NGC 3078 
is in fact a core galaxy. NGC 5576, on the other hand, has no detectable dust or strong radio 
emission. However, with 7maj = 0.36 and a = 0.79 it is close to an intermediate galaxy. NGC 3613 
and NGC 4365 are the only two core galaxies with disky isophotes: NGC 3613 harbours a thin 
edge-on nuclear disk inside the central arcsec, and NGC 4365 contains a kinematically decoupled 
core with the characteristics of a cold, rotating system (Surma 1993). 



5.4-2. Central cusp slopes 

In Figure we plot j'^^^^ as function of both the (major axis) break radius and the absolute 
B-hand magnitude Mb- Solid circles correspond to core galaxies, solid squares to power-law 
galaxies, and open circles to intermediate galaxies. All 10 core galaxies have break radii in the 
(relatively narrow) interval 0.3" ^ Tf, ^ 2", and are bright {Mb ^ —20). The power-law galaxies, 
on the other hand, typically have Mb ^ —20 and j'^^^ ^ 0.6. This once again confirms the 
dichotomy amongst early- type galaxies inferred from previous studies. 

But what about the intermediate galaxies? According to their isophotal structure, one would 



tend to classify them as core galaxies (see § 5.4.1). The same is true if one bases the classification 



on 7 (cf. Figure |6[). However, they also have properties that clearly differentiate them from 
core galaxies: their 7'-values are more similar to those of the power-law galaxies (although they 
occupy the low-7^g^j end of the distribution), their luminosity profiles have values of a that are 
significantly lower than the typical value for core galaxies, and finally, their average absolute 
magnitudes seem to be intermediate between those of both the core and power-law galaxie^: 
they can not be unambiguously classified as either one of those. More detailed investigations of 
their properties are required to see if they make up a truly, physically distinct class of objects, 
or whether these galaxies are distinct from both the core and power-law galaxies only because of 
distance and/or resolution effects. 



^except for MCG 11-14-25A, which is the faintest galaxy in our sample with AIb ~ 



-17.90 
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5.5. Misaligned Structures 

In fourteen galaxies we find radial intervals with strong isophote twists, non-zero fourth 
order moments (04 and/or 64) and local ellipticity maxima. A clear example is NGC 2699 
where at ~ 4" from the center the isophotes reveal an isophote twist of ~ 25°, both 04 and 64 
are positive, and there is a local maximum in ellipticity. In some cases the orientation angle 
of the misaligned structures changes with radius, giving it the appearance of spiral arms (e.g., 
ESQ 443-39, NGC 3595, and UGC 4551), wheras in other cases (ESQ 378-20, ESQ 443-39, 
ESO 447-30, NGC 2950 and UGC 6062) misaligned structures are evident in two distinct radial 
intervals. A clear example is ESO 447-30 with misaligned structures at both ~ 2.5" and ~ 10". In 
each of these galaxies, the two misaligned structures are not only misaligned with the host galaxy, 
but also with each other. 

Given the highly elongated nature of these structures and the fact that disks are common in 
elliptical galaxies, it is tempting to interpret these structures as stellar disks whose orientation 
axes are misaligned with those of the host galaxy. However, disks in spheroidal systems can only 
be stable if their symmetry axis is perpendicular to the equatorial plane of the potential, which 
does not result in a projected misalignment as observed. If, on the other hand, a stellar disk is 
embedded in a triaxial galaxy, projection effects can easily accommodate the observed misaligned 
morphologies. Furthermore, if the axis ratios of the triaxial host galaxy vary with radius, the 
projected angle of misalignment can vary with projected radius. Such configuration, therefore, can 
explain the multiple misaligned structures observed. A problem with this explanation, however, is 
that all 14 galaxies with misaligned structures are power-law galaxies which are generally thought 
to be spheroidal systems. Furthermore, if these structures are disks the observed fraction of disks 
embedded in power-law galaxies is too high (see § |5.6| below). 

A more likely explanation, therefore, is that these misaligned structures correspond to bars. 
This lends great support from the fact that isophotal structures, very similar to those seen here, 
have been observed in barred spirals and SOs (e.g., Wozniak et al. 1995; Friedli et al. 1996; Erwin 
& Sparke 1999). In particular, the bars within bars found by these authors provide a natural 
explanation for the multiple misaligned structures found in some of our galaxies. Indeed, both 
Kormendy (1981) and Wozniak et al. (1995) have previously noted the complex morphological 
structure of NGC 2950 (one of the fourteen galaxies in our sample with misaligned structures) 
and classified it as a double-barred galaxy. The bar interpretation is further supported by the 
fact that ten of the fourteen galaxies in our sample with misaligned structures are classified as SO 
or E/SO (see Table ||). It is unclear, however, how to interpret the misaligned structures in the 
four galaxies classified as ellipticals, although it is not inconceivable that these galaxies may be 
misclassified SOs seen close to face-on. Detailed investigations into the kinematics of these galaxies 
are required to confirm that indeed these structures correspond to tumbling triaxial potentials. 
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5.6. Stellar Disks 

In total 29 out of the 57 galaxies (51 percent) in our unperturbed sample reveal disky 
isophotes over at least a small radial interval (this includes the misaligned structures discussed in 
§ ^.5D . This diskiness is most straightforwardly interpreted as due to an embedded stellar disk 
(e.g., Scorza & Bender 1995, and references therein). 

The diversity of disk structures in our sample is large, ranging from extremely thin and 
small disks in the very center (i.e., NGC 3613) up to very large disks spanning the entire radial 
interval imaged (i.e., NGC 4621). In 18 galaxies the isophotes are disky in only a single radial 
interval, whereas the other 11 galaxies have disky isophotes in two radial intervals, separated by 
purely elliptical or boxy isophotes. This "double-diskiness" may be due to a nuclear disk inside 
the central hole of a separate outer disk or ring (i.e., van den Bosch &: Emsellem 1998), reflect the 
projected surface brightness contribution of a single disk which is locally dwarfed by that of the 
spheroidal, or be due to a bar-within-a-bar structure (see § ^.51) . The observed isophotal structure 
of an elliptical galaxy with an embedded stellar disk depends strongly on the ratio of the surface 
brightness distributions of the two components and on the inclination angle of the disk. Detailed 
disk-bulge decompositions (i.e., Scorza & Bender 1995; Scorza &: van den Bosch 1998) are required 
to decompose the observed surface brightness in the separate components, which is beyond the 
scope of this paper. 

Given the large fraction of early-type galaxies that reveal evidence for an embedded disk, and 
the fact that these disks are only visible if the disk is seen under a sufficiently large inclination 
angle, implies that the true fraction of early-type galaxies with embedded stellar disks is likely to 
be significantly larger than the 51 percent found here. In our unperturbed sample, 66 percent of 
all power-law galaxies have some amount of diskiness. This is a surprisingly high fraction. If we 
assume that all power-law galaxies have embedded disks, and that these disks are only detectable 
if their projected minor-to- major axis ratio is smaller than 0.5, one only expects to detect disks 
in 33 percent of all power-law galaxies (see Rix & White 1990 for a detailed discussion on the 
detectability of embedded stellar disks). The observed fraction of 66 percent implies either of 
the following: (i) all disks with a projected axis ratio smaller than 0.86 are detectable (assuming 
that all power-law galaxies harbor embedded disks), (ii) power-law galaxies are always seen under 
high inclination angles, (iii) the steep central cusp owes to the projection of the highly flattened 
structure (but see Lauer et al. 1995; Faber et al. 1997), or (iv) in some cases the diskiness is due 
to highly elongated prolate or triaxial structures (i.e., bars) rather than oblate disks. 

We only consider the latter of these explanations viable. It is supported by the fact that 
fourteen of our power law galaxies have their diskiness misaligned with the main galaxy. As we 
discussed in the previous section, these structures are most likely due to bars in SO galaxies. If 
we exclude these fourteen galaxies, we find embedded disks in only 32 percent of our power-law 
galaxies, which is in excellent agreement with the 33 percent expected if all power-law galaxies 
harbor embedded disks, but which are only detectable if the disk's projected axis ratio is smaller 
than 0.5. 
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5.7. Nucleated Galaxies 

Although the Nuker-law in general provides a good fit to the observed surface brightness 
of the central regions of early-type galaxies, in several cases the very centers reveal excess light, 
associated with an increase in slope as the HST resolution limit is approached. We follow Lauer et 
al. (1995) and refer to these features as nuclei. These unresolved nuclei may be due to non-thermal 
emission from AGNs or may reflect the presence of a nuclear star cluster. 

Table § lists the degree of nucleation observed, ranging from very mild, indicated by T (e.g., 
NGC 2634) to nuclei that stand out strongly against the 'background' Nuker-law profile, indicated 
by 'Iir (e.g., NGC 4482). In order to be termed 'nucleated' we require that the central steepening 
of the surface brightness occurs inside 0.15", and is visible from both the major and minor axes 
surface brightness profiles. Note that effects due to deconvolution as well as dust obscuration can 
produce artefacts in the luminosity profiles that are similar to our level I nucleation. In order to 
guard against such artificial 'detections', we have been fairly conservative in assigning nucleation. 

In total, nuclei are detected in 9 out of the 67 galaxies (13 percent) in our sample. This 
is less than half the detection rate of Lauer et al. (1995), who detected nuclei in 16 out of 45 
galaxies. Except for the fact that we have been fairly conservative in assigning nuclei, it is unclear 
why our detection rate is so much smaller. Only 2 of the 16 nucleated galaxies in the sample of 
Lauer et al. are core galaxies. In both cases the nucleation is interpreted as non-thermal emission 
from AGNs. In our sample, the two nucleated galaxies (NGC 4168 and NGC 5077) not classified 
as power-law galaxies have been detected at 1.4 GHz (see Table [l|). This is consistent with the 
idea that 'nucleation' in core (and intermediate) galaxies is associated with an active galactic 
nucleus. For none of the 7 nucleated power-law galaxies do we find any sign of non-thermal 
emission associated with an AGN (i.e., none of these galaxies have been detected at 1.4 GHz). 
This again is in good agreement with the results of Lauer et al. (1995). These results, however, 
are somewhat puzzling from an astrophysical point of view. An often advertised interpretation of 
(non-AGN related) nuclei is that they are stellar remnants of cannibalized, low-mass galaxies that 
have reached the center by means of dynamical friction. Since less luminous galaxies are in general 
denser than their more massive counterparts, one expects to find stellar nuclei predominantly in 
the less dense cores of the more massive galaxies, where tidal forces are less severe (e.g., Kormendy 
1984; Faber et al. 1997). The almost complete absence of stellar nuclei in core galaxies hints at 
some violent disruption process, most likely associated with massive black holes (see Merritt & 
Cruz 2001, and references therein). However, this fails to explain why power-law galaxies do seem 
to harbor stellar nuclei. These galaxies have massive black holes as well, and their intrinsically high 
stellar densities will only aid in disrupting any merger remnant that might reach the center of the 
potential well. A more viable explanation is that these nuclei are the product of some dissipational 
formation mechanism. The presence of highly fiattened disk-structures in the majority of these 
power-law galaxies supports the view that dissipational processes have played an important role. 
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6. Conclusions 

We have reported on the i2-band images of the central regions of 67 early-type galaxies 
obtained with the WFPC2 aboard the HST. This sample roughly doubles the number of early-type 
galaxies that have now been imaged at HST resolution, and complements similar data on the 
central regions of radio galaxies and the bulges of spiral galaxies. In this paper we have presented 
the data, focusing in particular on the isophotal structures and luminosity profiles. 

Our sample strikingly confirms the complex morphologies of the stellar component in the 
central regions of early-type galaxies which have become apparent from previous studies with 
HST. Whereas the outer regions of early-type galaxies can in general be classified as either disky 
or boxy, the isophotal structure in the inner ~ 10" varies strongly with radius, such that a bimodal 
classification does not suffice to describe the bewildering variety of structures seen. Instead, we 
have used a simple letter coding to describe the diverse morphologies. 

Previous imaging surveys with HST have shown that a more suitable criterion for classifying 
the central regions of early-type galaxies is based on their luminosity profiles (rather than their 
isophotal structure). Depending on the slope of the central cusp, galaxies are classified as either 
a core galaxy or a power-law galaxy. We determined the luminosity profiles along the major and 
minor axes for 57 galaxies in our sample (the unperturbed sample). The remaining 10 galaxies 
do not allow a proper assessment of their luminosity profiles because of dust obscuration. We 
have shown how luminosity profiles derived from an isophotal analysis of the data can produce 
relatively large systematic errors in the central region, owing to the sub-pixel interpolation used. 
We therefore use the actual pixel values in the central 0.2" to define the luminosity profiles. We fit 
the luminosity profiles with a Nuker-law profile (equation |Q]), which has been shown to adequately 
describe the luminosity profiles of early-type galaxies in the inner ~ 20". At radii less than 0.2" 
we fit this Nuker-law model directly to the observed pixel values, while at larger radii we use the 
luminosity profile as determined from the isophotal analysis. Detailed modeling has shown that 
this method allows an accurate recovery of the actual luminosity profiles. 

In a similar study, based on a sample of 61 early-type galaxies imaged with the WFPCl, 
Faber et al. (1997) classified galaxies as core galaxies if the central cusp slope 7 of the best fitting 
Nuker law was less than 0.3. Galaxies with 7 > 0.3 were classified as power-law galaxies. In our 
sample we find a subset of galaxies with luminosity profiles that have a continuously changing 
slope, rather than a double power-law shape with a clear break radius. For these systems 7 of 
the best-fit Nuker-law profile is smaller than 0.3, but it is only a measure of the gradient of the 
luminosity profile at radii much smaller than the resolution limit of the actual observations. At 
the resolution limit the actual slope of the luminosity profile is steeper than 0.3, and it is unclear 
whether to classify these systems as core or power-law galaxies. Although the isophotal properties 
are similar to those of the core galaxies, there is a hint that the luminosities of these galaxies are 
intermediate between those of the core and power-law galaxies. We therefore provisionally classify 
these galaxies as 'intermediate'. However, it is unclear at present whether these galaxies make up 
a physically distinct class or whether distance and/or resolution effects cause them to loose their 
distinct core and/or power-law characteristics. 
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F97 have pointed out that the distribution of central cusp slopes is clearly bimodal: whereas 
(by definition) core galaxies have 7 < 0.3, power-law galaxies are found to have 7 ^ 0.5. The 
region 0.3 < 7 < 0.5 is virtually devoid of galaxies. This bimodality closely follows the dichotomy 
in other properties: the core galaxies make up the bright end of the luminosity function and 
typically have boxy isophotes, whereas the power-law galaxies are fainter and typically have 
disky isophotes. We confirm this dichotomy with our sample: 93 percent of all galaxies that 
reveal disky isophotes are power-law galaxies, whereas 77 percent of all core galaxies have either 
boxy or purely elliptical isophotes. Furthermore, we confirm the bimodality of cusp slopes, but 
we do find a significant population of galaxies with 0.3 < 7 < 0.5. Combined with the class of 
intermediate galaxies, which have cusp slopes at the resolution limit that also fall in this range, 
we thus conclude that the bimodality is not as strong as previously claimed. 

Finally, the statistical properties of our sample can be summarized as follows: 

• In 15 percent of the galaxies in our sample we detect dust disks with diameters that 
range from 0.3" to 12.3". Together with the 28 percent of the galaxies in which we detect 
filamentary dust, the total detection rate of dust is 43 percent, in good agreement with 
previous estimates based on HST data (van Dokkum k, Pranx 1995). 

• In 9 galaxies we find evidence for nuclei. These are defined as (largely) unresolved, nuclear 

surface brightness spikes that cause a central upturn in the luminosity profiles. Our 
detection rate of nuclei (13 percent) is substantially less than the 36 percent found by Lauer 
et al. (1995) in their sample, but this is most likely a reflection of our more conservative 
approach to identifying nucleation. Whereas the nuclei in power-law galaxies are most likely 
stellar, the nucleation seen in core galaxies is generally non-thermal emission from an active 
galactic nucleus. 

• In 14 galaxies we find structures that are misaligned with the principal axes of the host 
galaxy. These structures arc most likely associated with stellar bars in SO galaxies. In 5 of 
these galaxies, multiple misaligned structures are apparent indicating bars within bars. 

• A large fraction (51 percent) of the early-type galaxies in our unperturbed sample reveal 
disky isophotes over at least a small radial interval, whereas this fraction increases to almost 
70 percent when only considering power-law galaxies. This is a remarkably large fraction: 
even if all power-law galaxies were to harbor an embedded stellar disk, one would not expect 
to detect these thin structures in more than ~ 30 percent of all cases because of inclination 
effects. If, however, we exclude the cases where the diskiness is misaligned with the main 
galaxy, and which are more likely associated with stellar bars, we find that ~ 32 percent of 
all power-law galaxies reveal some amount of diskiness. This is consistent with a picture in 
which all power-law galaxies harbor an embedded stellar disk. 

We are grateful to the referee. Tod R. Lauer, for his insightful comments that greatly 
improved the analysis and presentation of the data in this paper. The observations presented in 
this paper were obtained with the NASA/ESA Hubble Space Telescope, which is operated by 
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HF-01102.11-97.A awarded by the Space Telescope Science Institute. 
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A. Contour maps, luminosity profiles, and isophotal parameters 

In the following figures we present the data for all 67 galaxies in our sample. The upper 
right-hand panels show contour plots of the 33.0" x 33.0" field of view of the Planetary Camera. 
The thick black bar in the upper-right corner corresponds to an angular size of 5.0". Contours 
are plotted at intervals of 0.3 mag between fiji = 22.5 mag arcsec"^ and the surface brightness at 
r = 1.0" (in order to prevent central crowding of contours). 

The upper-panels on the left show the observed luminosity profiles (open circles are determined 
from the isophotal analysis, closed circles are the actual pixel values fitted and plotted as described 
in § [4.1| ) along the major (upper curves) and minor axis (lower curves, offset by +0.5 mag for 
clarity). The thick solid lines are the best fit Nuker-law profiles, the parameters of which are listed 
in Tables ^ and ^. In some cases data in the central regions is excluded from the Nuker-law profile 
fit because of nucleation and/or nuclear dust disks. In these cases a vertical thin dashed line 
indicates the lower limit of the data range used. 

The remaining panels plot the isophotal parameters for both the non-deconvolved image 
(solid triangles) and the deconvolved image (open circles with errorbars). For each galaxy we plot 
as function of radius the ellipticities (e), the position angles of the isophotes {9, measured from 
North to East), and the higher-order parameters and 6„ (with n = 3, 4). 

Some galaxies have their central regions saturated in the 700 sec exposure. In these cases 
we plot the isophotal parameters of the 300 sec exposures in the radial regime where the 700 
sec exposure shows saturation effects (rather than those of the combined image). In three cases 
(NGC 2950, NGC 3377, and NGC 4621), the central regions of the 300 sec exposures are also 
saturated. For these galaxies, only the isophotal parameters are plotted down to the radius where 
saturation effects in the 300 sec exposure start. 
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Fig. 1. — Isophotal parameters as a function of radius for a model galaxy which has perfectly 
elliptical isophotes with e = 0.26 independent of radius, and a central cusp with /(r) oc r~°''^ (cf. 
Figure Solid triangles correspond to the parameters of the (unconvolved) raw image, in which 
only Poisson noise is introduced. The open circles with error bars correspond to the parameters of 
the deconvolved image. Note that for r ^ 0.2" the parameters of the raw and deconvolved images 
are in excellent agreement. We therefore are confident that we can trust the isophotal parameters of 
our deconvolved images at least for semi-major axis radii larger than ~ 0.2". For r ^ 0.2", however, 
the sub-pixel interpolation of the fitting-routine systematically underestimates the ellipticity of the 
isophotes, and causes artificial deviatations from their perfectly elliptical shapes. 
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Fig. 2. — Characterization of errors in luminosity profiles. For each model i, the difference A//j(r) 
between the input surface brightness profile (i.e., equation [^) and the major axis luminosity profile 
determined from the isophotal analysis of the unconvolved image is calculated. For a sample of 8 
sets of model galaxies (each set is characterized by the same value of 7) the average (A//) (r) (right 
panel) and the corresponding standard deviation aAfj.{r) (left panel) is then computed, where the 
average is over the 50 models per set. The number in the legend corresponds to the characteristic 
central cusp slope 7 of the specific set. For r ^ 0.2" the luminosity profiles derived from the 
isophotal analysis are in excellent agreement with the intrinsic profile. For smaller radii, however, 
the error increases strongly, especially for the profiles with steep central cusp slopes. Note also 
that the average residual (A//)(r) strongly deviates from zero at small radii, implying that the 
main contribution to the error is of systematic rather than random nature. This is a reflection of 
the problems with the isophotal analysis due to the discrete, pixelized nature of the data and the 
sub-pixel interpolation. 
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Fig. 3. — In the upper panels, the sohd hnes correspond to the intrinsic major axis luminosity 
profile of the same model as in Figure [l|. Upper profiles (filled circles) correspond to the luminosity 
profiles determined from the isophotal analysis. Lower profiles (open symbols), offset by 1.5 magn, 
correspond to the luminosity profile determined with the new method utilizing the actual pixel 
values inside r < 0.2" (open squares) combined with the isophotal analysis data at r > 0.2" (open 
circles). The lower panels display the residuals for the two differents methods (same symbols as in 
upper panels). As is evident, the new method allows us to recover the intrinsic luminosity profile 
with high accuracy down to very small radii, whereas luminosity profiles determined solely from 
the isophotal analysis result in relatively large errors at small radii. 
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Fig. 4. — The left panel plots the average difference between the intrinsic cusp steepness 7modei 
and the best-fit Nuker-law profile 7fit for the set of models with 7 = 0.5 as a function of /. A 
value of / ~ 5 yields values for 7fit in good agreement with 7modei- For / ;^ 1, too much relative 
weight is given to the pixel values inside 0.2", and 7fit is too large. For too large /, is completely 
dominated by the isophotal luminosities outside 0.2" and the central cusp slope becomes poorly 
constrained, as is evident from the large errorbars. We therefore adhere to / = 5 for the final 
analysis. The right panel plots the fitted central cusp slope 7fit, determined with the new fitting 



method as described in § |4.1| , versus the intrinsic slope. The dots, which correspond to a single 
model, are grouped in 15 sets of 50 models each. Note that 7fit is in good agreement with 7modei- 
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Fig. 5. — Correlations between the two different measures of the central cusp of the luminosity 
profile: 7 is the asymptotic logarithmic slope of the best-fit Nuker-law profile, and 7' is the 
logarithmic slope of that fit at r = 0.1". The indices 'maj' and 'min' refer to the major and 
minor axes, respectively. Both 7 and 7' give very similar results along both axes as shown in the 
upper panels. Comparing 7 to 7' yields a more complicated pixture. Galaxies with 7 > 0.5 in 
general have 7' ~ 7. However, for 7 < 0.3 about half of the galaxies have 7' < 0.3, and the other 
half has 7' > 0.3. See the text for a more detailed discussion. 
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Fig. 6. — Histograms of 7 (left panels) and 7' (right panels) for the galaxies in our unperturbed 
sample. The indices 'maj' and 'min' refer to the major and minor axes, respectively. Black areas 
correspond to core galaxies, cross-hatched areas to power-law galaxies, and single-hatched areas 
to the intermediate galaxies. Although the distribution of cusp slopes reveals a clear hint for 
bimodality, we find less of a deficit of galaxies with 0.3 < 7 ^ 0.5 than F97. Note also that the 
intermediate galaxies themselves follow the same distribution of 7 as the core galaxies, but at the 
same time occupy the low-7' end of the power-law galaxies. 
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Fig. 7. — The parameter 7^^^^ as function of break radius ri, (left panel) and absolute i?-band 
magnitude Mb (right panel). Core galaxies are indicated by solid circles, intermediate galaxies by 
open circles, and power-law galaxies by solid squares. Note that the intermediate galaxies seem to 
have properties intermediate between those of core and power-law galaxies. 
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Table 1. Global parameters. 
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Table 1 — Continued 
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NGC 5173 


-19.52 


13.07 


32.8 


E/SO 


0.00 


1.043 


2419 




3.2 lb 0.6 


16.37 


NGC 5198 


-20.24 


12.42 


34.1 


E/SO 


0.14 


1.334 


2514 


0.845 


4.0 ±0.6 




NGC 5283 


-18.93 


14.13 


40.8 


SO 


0.10 


1.061 


3005 




13.4 it 0.4 




NGC 5308 


-20.13 


11.99 


26.6 


SOa 


0.82 


1.563 


1877 


0.739 






NGC 5370 


-19.45 


13.63 


41.3 


SO 


0.02 


1.145 


3053 








NGC 5557 


-21.52 


11.62 


42.5 


E 


0.21 


1.405 


3210 


0.824 






NGC 5576 


-20.00 


11.41 


19.1 


E 


0.23 


1.524 


1484 


0.807 






NGC 5796 


-20.38 


12.43 


36.5 


E 


0.25 


1.434 


2961 


0.953 


110 ±4 




NGC 5812 


-20.10 


11.86 


24.6 


E 


0.12 


1.420 


1961 


0.928 






NGC 5813 


-20.85 


11.11 


24.6 


E 


0.24 


1.663 


1917 


0.892 


16 ± 1 




NGC 5831 


-19.59 


12.06 


21.4 


E 


0.13 


1.397 


1660 


0.874 
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Table 1 — Continued 



Name 


Mb 




D 


Type 


£25 


d25 




B-V 


flA W-FIR 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) (11) 


NGC 5846 


-21.10 


10.74 


23.3 


E 


0.06 


1.631 


1722 


0.912 


22 ± 1 


NGC 5898 


-20.30 


11.80 


26.3 


E 


0.07 


1.459 


2171 


0.929 




NGC 5903 


-20.85 


11.58 


30.5 


E 


0.26 


1.530 


2512 


0.875 


32 ±2 


NGC 5982 


-21.32 


11.65 


39.3 


E 


0.30 


1.516 


2876 


0.816 




NGC 6278 


-19.70 


13.15 


37.1 


SO 


0.42 


1.345 


2791 






UGC 4551 


-18.77 


13.09 


23.6 


SO 


0.67 


1.322 


1722 


0.813 




UGC 4587 


-19.54 


13.49 


40.3 


so 


0.43 


1.215 


3060 






UGC 5467 


-18.87 


13.95 


36.7 


so 


0.00 


1.022 


2894 




7.9 ±1.1 13.40 


UGC 6062 


-18.82 


13.74 


32.7 


so 


0.39 


1.137 


2607 







Note. — Global parameters for the galaxies in our sample. Except for the data in the last two 
columns, these data are taken from the Lyon/Meudon Extragalactic Database (LEDA): 
Column (1): The name of the galaxy. 

Column (2): The absolute i?-band magnitude, based on the apparent magnitudes and distances 
listed in columns (3) and (4), respectively. 

Column (3): The total apparent i?-band magnitude corrected for galactic extinction and redshift 
effects (see Paturel et al. 1997). 

Column (4): Distances in Mpc, as derived from a pure Hubble expansion using the Virgo-centric 
infall corrected velocity (column [8]) and a Hubble constant of Hq = 80 kms~^ Mpc~^. 
Column (5): Morphological type, based on the code system of the Third Reference Catalogue of 
Bright Galaxies (de Vaucouleurs et al. 1991). A definition of this classification can be found in 
Paturel et al. (1997). 

Column (6): The ellipticity of the fiB = 25 mag arcsec"^ isophote. 

Column (7): The apparent diameter (in ^'^log of 0.1') of the galaxy according to the convention 
of the Second Reference Catalogue of Bright Galaxies (de Vaucouleurs et al. 1976). 
Column (8): The radial velocity (in kms^^) of the galaxy corrected for Virgo-centric infall (see 
Paturel et al. 1997 for the infall model). 

Column (9): Total B-V color, corrected for galactic extinction and redshift effects. 
Column (10): Radio flux (in mJy) at 1.4 GHz from the NVSS (Condon et al. 1998). 
Column (11): Far infrared magnitude, based on the IRAS 60 and 100 //m fluxes (see 
equation [|l|). 
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Table 2. Galaxy Morphology. 



Name 


Profile 


Dust 
Morph. 


Dust 
Level 


IC 


Nucl. 


M.S. 


(e) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


ESO 378-20 


\ 






dOd 




2 


0.34 ±0.02 


ESQ 437-15 




f 


III 


xd 






0.56 ± 0.08 


ESO 443-39 


\ 






OdOd 




2 


0.25 ± 0.04 


ESO 447-30 


\ 






OdOd 




2 


0.27 ± 0.08 


ESO 507-27 


\ 






dbd 


II 




0.49 ± 0.07 


ESO 580-26 


? 


f 


III 


X 






0.47 ±0.11 


IC 875 


\ 






do 






0.55 ± 0.05 


MCG 11-14-25A 


) 













0.11 ±0.02 


MCG 8-27-18 


\ 













0.10 ± 0.01 


NGC 2549 


\ 






dbd 


II 




0.41 ±0.11 


NGC 2592 


\ 


d 


0.6 


xOd 




1 


0.18 ±0.06 


NGC 2634 


\ 









I 




0.10 ±0.02 


NGC 2699 


\ 


d 


0.6 


xOdO 




1 


0.16 ±0.07 


NGC 2778 


\ 






OdO 






0.16 ± 0.05 


NGC 2824 


? 


f 


III 


X 






0.41 ± 0.20 


NGC 2872 


\ 


d 


0.8 


xO 






0.23 ± 0.01 


NGC 2902 


\ 


f 


II 


xO 






0.05 ±0.02 


NGC 2950 


\ 






dOd 




2 


0.29 ±0.04 


NGC 2986 


n 






?0 






0.18 ±0.01 


NGC 3065 


\ 


d 


0.4 


xOdO 






0.10 ± 0.06 


NGC 3078 


\ 


d 


1.2 


xb 






0.24 ± 0.01 


NGC 3193 


) 













0.20 ±0.03 


NGC 3226 


? 


f 


III 


X 






0.17 ±0.04 


NGC 3266 


\ 






OdO 


II 


1 


0.08 ±0.08 


NGC 3348 


n 













0.07 ±0.01 


NGC 3377 


\ 


f 


I 


d 






0.53 ± 0.01 


NGC 3414 


\ 


f 


I 


xO 






0.26 ± 0.02 


NGC 3595 


\ 






dbd 




1 


0.34 ± 0.04 


NGC 3613 


n 






dbO 






0.24 ± 0.06 
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Table 2 — Continued 



Name 


Profile 


Dust 
Morph. 


Dust 
Level 


IC 


Nucl. 


M.S. 


(e) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


NGC 3640 


) 













0.20 ± 0.03 


NGC 4121 


\ 






OdO 




1 


0.29 ± 0.04 


NGC 4125 


? 


f 


III 


xO 






0.34 ±0.12 


NGC 4128 


\ 






dbd 






0.39 ± 0.02 


NGC 4168 


n 


f 


II 


xO 


II 




0.18 ± 0.03 


NGC 4233 


? 


d 


12.3 


X 






0.30 ±0.18 


NGC 4291 


n 






Oh 






0.24 ±0.02 


NGC 4365 


n 






?dO 






0.29 ± 0.03 


NGC 4474 


\ 






dbd 


II 




0.20 ± 0.09 


NGC 4478 


\ 






db 






0.21 ± 0.07 


NGC 4482 


\ 






?0 


III 




0.27 ± 0.08 


NGC 4494 


? 


d 


1.6 


xO 






0.21 ±0.05 


NGC 4503 


\ 






OdO 




1 


0.26 ±0.01 


NGC 4564 


\ 






Od 






0.29 ±0.05 


NGC 4589 


) 


f 


II 


xb 






0.19 ± 0.03 


NGC 4621 


\ 






d 






0.37 ± 0.02 


NGC 4648 


\ 


d 


0.3 


OdO 




1 


0.21 ± 0.05 


NGC 5017 


\ 


d 


0.6 


xO 






0.17 ±0.02 


NGC 5077 


) 


f 


II 


xb 


II 




0.29 ±0.03 


NGC 5173 




f 


III 


x 






0.14 ±0.02 


NGC 5198 


) 













0.11 ±0.04 


NGC 5283 


? 


f 


III 


X 






0.13 ± 0.06 


NGC 5308 


\ 






dOd 






0.47 ±0.02 


NGC 5370 


\ 






Od 






0.18 ±0.11 


NGC 5557 


) 













0.20 ±0.03 


NGC 5576 


\ 






ObO 






0.29 ±0.03 


NGC 5796 


\ 


f 


I 









0.14 ±0.01 


NGC 5812 


\ 


d 


0.4 


xO 






0.02 ± 0.01 


NGC 5813 


n 


f 


II 


xO 






0.10 ±0.01 


NGC 5831 


\ 













0.30 ± 0.03 
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Table 2 — Continued 



NcLinG 


Profile 




UUST 


TP 


INUCl. iVi.o. 








Morph. 


Level 








(1) 


(2) 


(3) 


(4) 


(5) 


(6) (7) 


(8) 


NGC 5846 


? 


f 


III 


xO 




0.08 ± 0.03 


NGC 5898 


\ 


f 


I 







0.16 ±0.06 


NGC 5903 


n 


f 


I 


?b 




0.22 ± 0.03 


NGC 5982 


n 






Ob 




0.31 ±0.07 


NGC 6278 


\ 






bOd 


1 


0.22 ± 0.09 


UGC 4551 


\ 






Od 


I 1 


0.15 ±0.04 


UGC 4587 


\ 


f 


I 







0.28 ±0.07 


UGC 6062 


\ 






OdOd 


2 


0.30 ± 0.04 



Note. — Isophotal parameters for all the 67 galaxies in our sample. Column 
(2) indicates the classification where \, ), and fl indicate power-law, intermediate 
and core galaxies. A question mark indicates that no meaningful luminosity 
profile could be determined (either due to dust or nucleation) and the galaxy 
classification is therefore unknown in these cases. Column (3) gives the dust 
morphology where 'd' stands for dust disk and 'f stands for dust filaments. 
Column (4) indicates either the amount of filamentary dust ranging from I (small) 
to III (large) or the diameter of the dust disk in arcsec, if applicable. Column 
(5) gives the isophotal code (IC) of the galaxy going from r = 0.2" to the outer 
region where d, b, 0, x, and ? indicate disky, boxy, neither disky nor boxy, 
undetermined due to dust, and undetermined due to small surface brightness 
gradient. Column (6) indicates the degree of nucleation, ranging from I (weak) 
to III (strong). Column (7) indicates if and how many misaligned structures 
(M.S.) are detected in the galaxy (see § 5.5). Column (8), finally, lists the mean 
ellipticity (determined over the radial interval 1.0" < r < 10.0") and its standard 
deviation. 
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Table 3. Parameters of major axis luminosity profiles. 



Name 


Profile 


h 
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TTcin Pisn 9fi 


? 














Tr" S7c; 

ivj o / 


\ 
\ 


1 7 fit; 

1 ( .DO 


U.oO 


1 81 
1.01 


1 (17 
l.U 1 


1 87 
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1 f! R1 
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IVPP 9fiQQ 
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Z.DD 


n 81^ 

U.oO 
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1 7 na 
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n Q/i 
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n QQ 
u.yy 
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INvjv^ ZoZ4 
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\ 


18 11 
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l.Ul 
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\ 


1D.041: 


1 

i.yo 


1 ^^7 
1.0 ( 


n /1Q 


n 88 

U.OO 


n f^n 

U.OU 


NGC 2950 


\ 


16.16 


2.40 


1.81 


0.82 


2.43 


0.82 


NGC 2986 


n 


16.09 


1.77 


1.50 


0.18 


1.24 


0.20 


NGC 3065 


\ 


17.97 


0.99 


2.28 


0.79 


3.14 


0.84 


NGC 3078 


\ 


17.19 


2.35 


1.60 


0.95 


3.43 


0.95 


NGC 3193 


) 


15.46 


0.59 


1.89 


0.01 


0.81 


0.44 


NGC 3226 
















NGC 3266 


\ 


17.78 


1.34 


2.06 


0.64 


2.25 


0.66 


NGC 3348 


n 


16.00 


1.17 


1.53 


0.09 


0.99 


0.18 


NGC 3377 


\ 


15.37 


3.15 


1.35 


1.12 


1.42 


1.12 


NGC 3414 


\ 


16.46 


1.40 


1.45 


0.83 


1.72 


0.84 


NGC 3595 


\ 


17.43 


2.33 


1.52 


0.75 


2.30 


0.76 


NGC 3613 


n 


15.11 


1.53 


1.06 


0.04 


0.34 


0.17 
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Table 3 — Continued 



Name 


Profile 


h 


a 




7 
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7' 
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(2) 
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(5) 
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z.Uz 


n 1 n 

u.iy 


INvjU 4Zoo 


7 














ivir^r^ /I0Q1 
iNVjU 4zyi 


1 1 


10. UD 


1 Qc; 

i.OO 


1 fiO 
l.DZ 


n nn 
U.UU 


n An 
U.OU 


n 1 Q 
U.io 


IN^jU 4oDO 


1 1 


1 HQ 

iO.Uo 


OQ 

z.zy 


1 OQ 
l.ZO 


n 1 Q 
U.io 


1 CO 
l.OZ 


n 1 Q 
U.io 


l\Tn AA7A 


\ 


i ( .y ( 


O OQ 


1 RK 
i.OO 


n 7o 
U. ( z 


A KK 
4.00 


n 7o 
U. / z 


ivrr^ r^AA vs 


\ 


i / .0 / 


u.yo 


1 7Q 
i. / O 


n p;7 
U.o / 


/I 1 A 
4.10 


n Ai 
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U.Oi 


Air^r^ AKQCi 
L\(jVj 4ooy 


) 


1 RA 


n Qo 


1 /I /I 
i.44 


n o/^ 
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i.iZ 


1 /I Q 

1.4o 


1 1 o 
I.IZ 




) 


1 R /IQ 


1 no 
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n OQ 
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i.Oi 


n Qn 
U.oU 


ivpp p;i7^ 

IN 01 / O 


7 














NGC 5198 


) 


14.79 


2.61 


1.13 


0.23 


0.16 


0.42 


NGC 5283 


? 














NGC 5308 


\ 


15.55 


0.39 


1.27 


0.82 


0.83 


0.96 


NGC 5370 


\ 


17.84 


0.78 


1.50 


0.62 


1.63 


0.71 


NGC 5557 


) 


16.18 


0.80 


1.77 


0.14 


1.21 


0.33 


NGC 5576 


\ 


15.79 


0.79 


1.73 


0.36 


1.41 


0.51 


NGC 5796 


\ 


16.19 


0.83 


1.67 


0.41 


1.06 


0.56 


NGC 5812 


\ 


16.61 


1.25 


1.67 


0.59 


1.84 


0.62 


NGC 5813 


n 


15.77 


1.90 


1.33 


0.00 


0.73 


0.03 


NGC 5831 


\ 


16.90 


0.47 


1.84 


0.33 


1.78 


0.64 
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Table 3 — Continued 



Name 


Profile 


h 


a 


/3 


7 


n 


7' 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


NGC 5846 


7 














NGC 5898 


\ 


16.32 


1.23 


1.57 


0.41 


1.24 


0.46 


NGC 5903 


n 


16.84 


1.80 


1.48 


0.13 


1.59 


0.14 


NGC 5982 


n 


15.60 


0.99 


1.52 


0.00 


0.74 


0.18 


NGC 6278 


\ 


15.69 


0.76 


1.62 


0.55 


0.60 


0.77 


UGC 4551 


\ 


17.06 


2.19 


2.16 


0.51 


2.26 


0.51 


UGC 4587 


\ 


15.48 


1.04 


1.24 


0.72 


0.19 


0.90 


UGC 6062 


\ 


18.35 


0.90 


1.81 


0.80 


2.75 


0.85 



Note. — Parameterizations of the major axis luminosity profiles. Columns (2) indicates the 
classification of the galaxy (copied from column (2) in Table ^). Columns (3) - (7) list the best fit 
parameters of the Nuker-law (equation |^). Here the break radius, r^, is in arcsec and in mag 
arcsec"^. Column (8) lists the value of 7' (see § ^ for definition). 
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Table 4. Parameters of minor axis luminosity profiles. 



Name 


Profile 


h 


a 




7 


n 


y 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


H/Ow oio-ZU 


\ 


I 7 Qfi 

I I .yo 


u.oy 


9 nfi 


n 87 

U.o I 


1 Q8 
i.yo 


11c; 
1.10 




? 
















\ 
\ 


1 7 ^9 
i / .OZ 


1 .00 


i.DO 


n 71^ 

U. ( 


1 1 Q 
i.iy 


n 77 


"pcin /I/17 "^n 


\ 
\ 


1 P. 88 
iO.OO 


1 '^1 
i.Oi 


1 fin 

i.DU 


n 8"^ 

U.oO 


1 1 Q 
i.iy 


n 8fi 

U.oD 


"PQn 1^07 97 


\ 
\ 


1 7 71 


i.OO 


9 91^ 

z.zo 


n fi7 


1 fifi 

l.DD 


n fiQ 
u.oy 


TTcin Pisn 9fi 


? 














Tr" S7c; 

ivj o / 


\ 
\ 


io.yo 




1 78 
1. 1 o 


1 H/l 


9 1/1 


1 91 
l.Zl 


A/TPP 11 1/1 9^ A 
lViV_/V^ li-14-ZO-rv 


^ 
) 


io.oy 


n 7n 

U. 1 u 


9 91 
Z.Zi 


n nn 
u.uu 


n ^fi 

U.OO 


n ^1 

U.Ol 


A/TPP 8 97 1 S 


\ 
\ 


1 8 "^Q 

io.oy 


n 89 

U.oZ 


1 87 
i.o / 


n 7Q 
u. 1 y 


1.00 


n 80 
u.oy 


MPP 91^40 


\ 


1 ^ 88 
10. oo 


o.4ty 


1.00 


U.OO 


n QzL 


U.OO 


INvjr^^ ZOyZ 


\ 
\ 


1 7 /in 


9 /in 


1 8fi 
i.oD 


1 nfi 


1 09 

i.yz 


1 nfi 

l.UO 


1\TPP 9f?Q/l 


\ 


1 7 

i t .OO 


Z.OD 


1 K/l 


n 70 
u. < y 


1 fif? 

l.DD 


n 8n 

U.oU 


IVPP 9fiQQ 


\ 


1 7 87 


i.UD 


1 8/1 


n Qc; 
u.yo 


9 /in 


n 08 
u.yo 


TvTPP 9778 
IN vjtl^ Zt to 


\ 


ID. 00 


1 /IQ 


1 K'i 
l.Oo 


n Of? 
u.yo 


n 77 


n OB 
U.yo 


MPP 989/1 
INvjv^ ZoZ4 


7 














TMPP 9879 
IN vjr^ Zo 1 Z 


\ 


17 8/1 


^19 
O.IZ 


1 1^8 
l.Oo 


1 nn 

l.UU 


9 7fi 
Z. ( 


1 nn 

l.UU 


MPP 9Qn9 
INvjrV^ ZyUZ 


\ 


1 97 


9 1/1 


1 t^n 

l.OU 


n /1Q 


n 77 


n f^i 

U.Ol 


NGC 2950 


\ 


16.16 


1.66 


1.78 


0.84 


1.65 


0.84 


NGC 2986 


n 


16.05 


1.86 


1.45 


0.19 


0.97 


0.20 


NGC 3065 


\ 


18.04 


0.86 


2.17 


0.83 


2.97 


0.89 


NGC 3078 


\ 


17.49 


3.42 


1.66 


1.06 


3.26 


1.06 


NGC 3193 


) 


14.21 


1.18 


1.18 


0.00 


0.14 


0.48 


NGC 3226 
















NGC 3266 


\ 


17.28 


2.39 


1.91 


0.67 


1.46 


0.67 


NGC 3348 


n 


16.06 


1.08 


1.61 


0.06 


0.98 


0.18 


NGC 3377 


\ 


15.91 




1.29 








NGC 3414 


\ 


16.29 


1.13 


1.39 


0.85 


1.10 


0.89 


NGC 3595 


\ 


17.26 


2.43 


1.75 


0.77 


1.34 


0.77 


NGC 3613 


n 


15.73 


1.23 


1.46 


0.13 


0.65 


0.25 
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Table 4 — Continued 



Name 


Profile 


h 


a 


P 


7 


n 


y 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


OU4U 


) 


1 c; /1 1 

10.41 


U.OD 


1 OS 
l.ZO 


n 1 7 

U.l ( 


n /IS 

U.40 


n /in 

U.4U 




\ 
\ 


1 n Qo 
Ly.oz 


i.zo 


o nn 


n QK 
U.OO 


o.4o 


n 

U.o( 


InvjU 41Z0 


7 
















\ 


i ( .oU 


U.Oo 


00 

z.zz 


n A7 
U.D ( 


1 OS 

i.yo 


n 70 
U. / y 


InvjU 4iDo 


1 1 


i / .40 


i.4o 


1 Qn 
i.oU 


n 01 


1 A/i 

1.D4 


n OQ 
U.Zo 


INvjU 4Zoo 


7 














ivir^r^ /I0Q1 
iNVjU 4zyi 


1 1 


10. U ( 


i.oD 


i.DO 


n nn 
U.UU 


n /i7 
U.4/ 


n 1 s 
U.io 


IN^jU 4oDO 


1 1 


iO.ZO 


1 KA 
1.04 


1 QO 

i.oy 


nil 
U.ii 


1 t;c; 
i.oo 


n 1 Q 
U.io 


l\Tn AA7A 


\ 


1 V /I Q 

i ( .4o 


Q QQ 
O.OO 


i.yo 


U.OO 


QS 

z.oo 


n KQ 

U.OO 


ivrr^ r^AA vs 


\ 


i ( .Do 


U.o / 


i.oi 


n t;/i 
U.o4 


Q Q1 

o.oi 


n c;o 

u.oy 


INLiLy 44oz 


\ 
\ 


1 o no 


O 1 V 


1 KCi 
i.OU 


U.OU 


A no 
4.Uz 


n c;n 
U.OU 


ivr^r^ AAQA 

INvjU 44y4 


7 


io.oo 




1 1 Q 

i.io 








IN^jU 40Uo 


\ 


io.oy 


1 01 

l.zl 


1 An 
i.OU 


n AQ 

U.OO 


1 7n 
i. ( U 


n AA 
U.OD 


rs\ji\j 4004 


\ 


1 on 
io.yu 


1 1 Q 

i.io 


1 OQ 

i.yo 


n 7Q 
U. / 


1 n 
z.iu 


n so 
U.OZ 


i^KjLj 4ooy 


) 


1 K RA 
10.04 


U. ( i 


1 /It; 
1.40 


n 1 Q 
U.io 


n KA 
U.o4 


n /I7 
U.4 ( 


'KJCr^ A ROT 


\ 


1 V KO 

i / .oz 


U.Dy 


1 /in 
i.4y 


u.yo 


/Id 
4.oi 


1 nn 
i.UU 


INljLy 4D4o 


\ 


io.oU 


o.iu 


1 yt; 

i. ( 


1 HQ 

i.Uo 


1 on 
i.zy 


1 nQ 
i.Uo 


IN IjU OUi ( 


\ 


11.10 


1 Q1 

i.oi 


i.UO 


1 1 n 

i.iy 


O 1 o 
Z.lZ 


1 1 o 

i.iy 




) 


ID. 40 


lie; 
i.iO 


1 7/1 
i. (4 


n oc; 
U.zo 


1 1 o 

i.iy 


n QQ 
U.OO 


ivpp p;i7^ 

IN 01 / O 


7 














NGC 5198 


) 


14.81 


2.04 


1.15 


0.21 


0.14 


0.52 


NGC 5283 


? 














NGC 5308 


\ 


18.39 


0.85 


2.18 


0.86 


3.95 


0.92 


NGC 5370 


\ 


18.61 


1.26 


2.51 


0.64 


2.40 


0.68 


NGC 5557 


) 


15.74 


1.11 


1.54 


0.18 


0.67 


0.33 


NGC 5576 


\ 


15.57 


0.91 


1.76 


0.36 


0.89 


0.53 


NGC 5796 


\ 


16.35 


0.82 


1.71 


0.40 


1.05 


0.57 


NGC 5812 


\ 


16.69 


1.22 


1.73 


0.58 


1.91 


0.62 


NGC 5813 


n 


15.88 


1.61 


1.41 


0.00 


0.74 


0.05 


NGC 5831 


\ 


15.16 


0.99 


1.34 


0.39 


0.27 


0.65 
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Table 4 — Continued 



Name 


Profile 


h 


a 




7 




7' 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


NGC 5846 


? 














NGC 5898 


\ 


16.63 


0.90 


1.61 


0.38 


1.35 


0.49 


NGC 5903 


n 


16.66 


2.30 


1.35 


0.12 


0.98 


0.13 


NGC 5982 


n 


15.26 


1.91 


1.42 


0.05 


0.46 


0.12 


NGC 6278 


\ 


16.31 


1.11 


1.96 


0.65 


0.86 


0.76 


UGC 4551 


\ 


17.39 


1.95 


2.57 


0.49 


2.48 


0.49 


UGC 4587 


\ 


15.97 


0.72 


1.47 


0.63 


0.23 


0.92 


UGC 6062 


\ 


18.46 


1.05 


2.05 


0.80 


2.16 


0.85 



Note. — Same as Table ^, except now for the minor axis luminosity profiles. 



Table 5. 



Profile 


d 


db 


b 





Total 


Power-law 


20 


7 


2 


12 


41 


Intermediate 








2 


5 


7 


Core 


1 


1 


3 


4 


9 


Total 


21 


8 


7 


21 


57 



Note. — Correlation of luminosity profiles with isophotal structure. For each class of luminosity 
profiles (power-law, intermediate, and core) we indicate the number of galaxies (in the unperturbed 
sample) with a common isophotal characteristic. We distinguish between disky 'd', boxy 'b', disky- 
boxy 'db', and regular '0' (see § ^.4.1 for definitions). Note that most disky and disky-boxy galaxies 
are power-law galaxies, whereas most boxy galaxies are classified as either core or intermediate 
galaxies 



